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Abstract 

A model Hamiltonian is proposed in order to understand the localization-delocalization 
transition in a quantum dot, where there are two gate voltages: top and side. Considering 
energetically favorable degrees of freedom only, we achieve a finite dimensional Hilbert space. 
As a result, exact diagonalization is performed to find the ground state energy of the system. 
It is the purpose to explain the peculiar pattern of the electron addition energy measured 
in the dot of two gate voltages. 
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Many-body effects in low dimensional systems have attracted many interests in both 
experimental and theoretical point of view. As one of low dimensional systems, quantum dots 
are fabricated and used in order to investigate notorious problems of electron correlation. 
Ashoori invented the useful tool called single electron capacitance spectroscopy|Jl|] to measure 
electronic properties of quantum dots. It became possible to measure directly the iV-electron 
ground state energies of quantum levels of a dot as a function of magnetic field0. Further 
studies on electron addition spectra of quantum dots showed that there are bunches in 
electron additions!^. This strange electron correlation of bunching was investigated by 
intensive theoretical efforts^, ||, |^. 

Recently, Ashoori group observed the localization-delocalization transition^, using newly 
designed dots with two kinds of gate voltages: top Vj^ and side Vs. The side gate voltage 
plays a crucial role in analyzing the edge state localization. In this seminal experiment, what 
they measured is the dependence of the capacitance peak on and Vs in electron additions. 
They plotted the lines of versus Vs for each A^-electron in the dot up to about = 100. 
The features of the plot are summarized as the followings. a) Two kinds of lines are observed: 
one has a small slope, and the other is steep. b)For the low densities of electrons, the spacings 
between adjacent lines are irregular, c) However, one observes the general trend of decreasing 
of the spacings as more electrons are added. d)The small-slope lines gradually become more 
steep as the number of electrons in the dot is increased. e)The anticrossing takes place when 
a small-slope line and a steep line are merging. The interesting observation connected with 
the anticrossing is that there are two kinds of anticrossings: normal and abnormal. The 
anomalous anticrossing shows that the chemical potential of A^ + 1-electron state is lower 
than that of A^-electron state. In fact, this striking result shows that the edge localized 
electrons appear to bind with electrons in the dot center. 

It is the purpose of this Letter to explain the localization-delocalization transition with a 
tractable model Hamiltonian. A reduction of the corresponding Hilbert space is proposed in 
order for the system to be calculable. This truncation is called single level approximation, 
which is resemblance of the lowest Landau level approximation in the fractional quantum 
Hall effect ip, ^. Up to the single level approximation, we diagonalize the Hamiltonian 
exactly, using the Lanczos method. It is found that the spacings between small-slope lines 
are attributed by Coulomb blockade. It is understood that the steep lines are related with 
localized states. Furthermore, we notice that the anomalous anticrossings are possible by 
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the quantum interference. In consequence, all features of the plot described in the above, 
except for b), are explained in this Letter. 

In order to study the quantum dot, which is experimentally investigated in Ref. 7, we 
consider a model Hamiltonian, which is written as a function of the side gate voltage Vs: 

(1) m)=^ext + ^loc + ^int' 

where -ffext describes the inside of the dot relating with extended states, H^q^, corresponds to 
localized states, and H^^^ represents interactions between extended electrons and localized 
electrons. The measured top gate voltage will be a function of Vs in relation with the 
electron addition energies: aeAVl{N) = Eg{N + 1) - 2Eg{N) + Eg{N - 1), where ^g(A^) 
is the ground state energy of H{Vs) for total iV-electron in the system, and the parameter 
a is a geometrical coefficient. 

While it is difficult to present H-^q^ in terms of position and momentum variables, we 
write the extended state Hamiltonian -ffext -B-field as 
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(2) ^ext = E{AlP"'^ + -^^l' + 9"^*^'l/^"'^l'-^^^-^}+ E 

^-^ 2m* c 2 m*c ^-^ k p,- — 

where we adopt a two-dimensional pancake-type dot with a confining potential controlled 
by uj, and the position of the i-th electron is presented by pi. The (^-factor is usu- 
ally given by a very small value enough to ignore the spin term. The dielectric con- 
stant K is introduced in the Coulomb interaction. The free part of -f^gxt solved by 



Fock|]10[. In fact, introducing one-particle creation operators cj^^^, we find the corresponding 



eigenfunctions: ^„„^(p,0) = < p\cl^JO >= ^2A/'K{{n - 1 - |m|)/2)!/((n - 1 + |m|)/2)! 
exp(— zm</)) exp(— y4p^)(2Ap^)l™'l/^ L^^^_-^^_^^^^y^{2Ap'^) x^^^^ where the principal quantum 
number n runs as 1, 2, 3, ■ ■ ■; the corresponding magnetic quantum number m is given 
by — (n — 1), — — 3), ■ ■ ■ , (ra — 3), {n — 1); the spin index a =t, i; and the short 
hand notation A = m*u/2h, oj = ^cj^ + \uJq, luq = eB/m*c. Here, following the normal- 
ization of Arfken|jTl|, we note the associated Laguerre polynomials: L'^{x) =J2m=o(~^)"^ 
{n + k)\/{{n — m)\{k + m)\m\}x"^ = X]m=o ^n("^)-^"^' where the fractional numbers of y4^(m) 
will be used to present the coefficients for the Coulomb interaction in the formalism of second 
quantization. 

It is a straightforward process to obtain the second quantized Hamiltonian. With com- 
plicated coefficients for the Coulomb interaction, the second quantized Hamiltonian will be 
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written as the usual form in terms of c}^^^ and Cnma- For this full Hamiltonian, the dimension 
of the corresponding Hilbert space is infinite. It is impossible to calculate an exact ground 
state energy. Thus, we need truncation. In the case of a small S-field, the principal quan- 
tum number n plays the role of distinguishing shells. When we consider electrons in the 
dot, we note that the electrons occupy from lower energy states. There will be the biggest 
value of n in this situation. We can divide the number of electrons as N — A^core + -^shell' 
where A^core — {n — l)n, and -/Vgjig}} = 0, 1, • • • , 2n. In this Letter, reducing the degrees of 
freedom, we ignore detailed interactions between principal quantum levels, and also neglect 
higher energy states. This is called single level approximation. Like the case of Ref. 7, we 
simply let the B-ficld zero from now on, hence Uc = and u = u. Roughly taking care of 
the interaction between the core and the shell electrons, we introduce parameters C(iVcore)- 
In consequence, the truncated extended state Hamiltonian is written as 

(3) i/ext = ^(^core) + C^(iVcore)iVop + hunNop 

— {n-l)<k,l+m,l,k+m<{n—l) cr,a' 

where A'op = X^mo- ^nma^nmai ^ud the zero-point energy E{Ncore) effectively describes the 
interaction among the core electrons. Now, the dimension of the corresponding Hilbert 
space is finite, in fact, the number of different ways in taking A^gj^ell out of 2n. Thus, it 
is calculable. The values of E{Ncore) and C(A^core) will be determined later. Using the 
simple formula: l/|p2 - pi\ = J2meZ exp(?m(0i - 02)) dkJm{kpi)Jm{kp2), where Jm{x) 
is the Bessel function, we calculate the coefficients Vn{k, Z, m) of the Coulomb interaction in 
the level of n: 



(4) K(A:,Z,m) = — VI 



n—l—\l+m\ I n— 1 — |i| | n—l — \k+m\ | n— 1— |fc| | 
2 2 • 2 • 2 • 

[ n-l+|/| [ n-l+\k+m\ [ n-l+|fc| [ 
2 ■ 2 ■ 2 ■ 2 ■ 



^ ^ L-l-|i+m| (^l)^ n-l-lil (fe2)^ L-l-jfc+m| (^3)^ Lll-|);| (^4) 



w-l-|i+m| w-l-|i| n-l-\k + m\ n-l-|fc| 
2 2 2 2 



ki=0 k2=0 k3=0 k4=0 



x{ki + k2 + ^)K^3 + k4+ ^)! 



fci+fc,+ l'+'"l+l'l-KI fc3+fc^+ lfe+'"l+lfc|-l^ 



h=0 l2=0 ^ ^ 

r(|m| + /i + /2 + 5) 



2\m\+h+l2 
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Note that the overall factor is proportional to ^ as e^yfAjlK = V — Voy^hu/lmeV. 
Comparing y/cu with hcu, which is the energy difference between principal quantum levels, 
we find that our single level approximation is the more valid for the larger value of cu, that 
is, the stronger confinement. 

Turning our attention to localized electrons, we consider only the case of a single localized 
state without loss of generality. We guess the Hamiltonian H^^^ with the creation operator 
dl for the locahzed electron as 

(5) Hi^^ = e(n^ + + Un^rii, 

where = dj^d^. It seems that two localized electrons at the same site are most likely feel 
a large effect of repulsion. Thus, the value of U should be large. The energy value of 2e + C/ 
is energetically unfavorable. Thus, it is enough to consider only two cases, or e, for the 
energy of if^Q^,. 

The interaction Hamiltonian between extended and localized electrons is written as 

(6) -^int ^ ^ {^nm.clima^nmai'^l ~l~ ^J.) ~l~ ^nm^nma'^T ~l~ ^nm^a^nma} ■ 

n,m,(7 

Here we consider the direct Coulomb interaction and the tunneling effect. Since the localized 
electron wave function < p|(ij.|0 > is not known, we can not calculate A„m, nor tnm- 

We have introduced the Hamiltonian of the system, H[Vs). Our task is now to find the 
ground state energy of H{Vs). The strategy for this is to calculate the ground state energy 
of i^ext + -^loc fi^st) Rayleigh-Schrodinger perturbation theory with respect 

to i^in^. And then, we make connection between and Vs, using the electron addition 
energy. 

We have already determined the energy of H^q^ trivially. We consider i?ext ^ corre- 
sponding principal quantum number n with the single level approximation. We write the 
ground state energy of i^ext (3) ^ 

(7) E{N) = E{Ncore) + C{Ncore)N^i,e\l + ^^^^shell + ^c(iV), 

where Ec{N) is the Coulomb correlation energy, which plays an essential role in the electron 
addition energy. In the single level approximation, it is obvious to notice i?c(^core + 1) = 0. 
We calculate the correlation energy Ec{N) up to A'" = 170 by using the Lanczos method, for 
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instance, Edl) = 0.0, Eci2) = 1.772451/, Ec{S) = 0.0, Ec{A) = 0.88622^, • • • , Ec{156) = 
129.724^, • • • , and £;c(170) = 39.8113^. 

Using the ground state energy of -ffext (^)' calculate the electron addition 

energy E{N + 1) - 2E{N) + E{N - 1) as 



(8) 



h(jj + A{n) for = n(n + 1) , 

Ec{n{n + l) + 2) for AT = n(n + 1) + 1, 

£;c(A^ + 1) - 2EciN) + EciN - 1) otherwise. 



where A(n) = C{n{n + 1)) - C((n - l)n) - Ec{n{n + 1)) + Ec{n{n + 1) - 1). Since it is 
numerically shown that the value of Ec{n{n + 1) — 1) — Ec{n{n + 1)) contains a negative 
factor proportional to n, the parameter C{n{n + 1)) must cancel the factor by subtraction 
so that hw + A(n) is positive for all n to follow the concept of Coulomb blockade. The 
parameters C{n{n+ 1)) are determined by A(n), which should be properly chosen to satisfy 
experimental data. All values of E[Nqqyq) and C(A^core) are determined recursively from 
E(0) = C(0) = 0. 

We write two candidates I^'at+i > and d'll^ ^ > for the ground state of the Hamiltonian 
ifext+-f^loc -^+1 electrons as (ifext+-f^loc)l*^+i >^ -B(A^+l)|*iv+i >, and (i?ext + 
-^loc)^^-!^^ ^)'^V\^ N >■ Considering the total Hamiltonian of Eq. (1) now, we 

calculate the correction energy of the first order perturbation as < \l/Ar+i|i/jj^^|^jv+i >= 

and < ^iv|rf^i7jj^^4l^Af >= < N\Y.n,m,a^nmCimaCnma\^ N > = A(A^). NotC that the 

ground state energy of the system is E{N + 1) or E{N) + e + A(iV). If these two values 
are almost same, then this is the degenerate case and we should diagonalize a matrix in 
perturbation. With the two degenerate states d^l^Ar > and I^at+i >, we find the relevant 
Hamiltonian i^cross^ 



(9) 



cross 



E{N)+e + X{N) t*{N) 
t(N) E(N+1) 



where < N+i\^nm a'^nrn<^yima\'^ n >= . We find that the degeneracy is removed 
by t{N), and the ground state energy is given by 1{E{N + 1) + E{N) + e + X{N) 
—s/(E(N + 1) - E{N) - e - X(N)y + 4|i(A^)|2} = £;cross(A^ + 1)- In consequence, the 
ground state energy of the system Eg{N + 1) up to the first order perturbation with 
respect to H^^^ is given by min{£'(7V) + e + X{N),E{N + 1)} or Ectoss{N + 1) if 
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E{N) + e + A(A^) ^ E{N + 1). This result of Eg{N + 1) and also Eg{N) will be used 
in the below. 

Since a line measured in the experiment presents the event of single electron oscillation 
between the dot and the contact, it is appropriate to use the notation of ^ in eV^{N + |). 
Introducing the offset value e\\{^), which is the first line in the plot of \\ versus Vs, we find 
eVtiN + i) = eVt(|)+ eAVt(z) = eVt(i)+ ^{Eg{N + 1) - Eg{N)- Eg{l) + Eg{0)}. 
Note that the energy differences D(N + 1) = E(N + 1) - E(N) - X(N) satisfy the inequality 
of D(N) < D(N + 1) for all N with small A(A^). We find that Eg(N + 1) - Eg{N) in 
eVi{N + I) is given by one of the five expressions according to e: 

D(N) + X(N) ioTe<D(N), 

E{N) + e + X{N) - Ecross{N) for e ^ D{N), 

(10) {e + X(N) for D(N) < e< D(N + 1), 
EcrossiN + 1) - E{N) for e !^ D{N + 1), 
D{N+1) + X{N) for e> D{N +1). 

As far as parameters in the Hamiltonian are functions of Vs, Eg{N) is also a function of Vs- 
Hence, we have connected Vg to V^{N + ^). In Rcf. 7, observing the large capacitance shows 
that the localized states exist at the periphery of the dot. The energy value of the localized 
state e is much affected by Vs. Furthermore, the side gate voltage Vs effectively changes the 
dot confining potential. Thus, we assume the dependence of the side gate voltage on the 
parameters as 

(11) e^eo + Pe\Vs\, hu; ^ huo + -fe\Vs\, eVt(^) = eVto + ^el^sl- 

Summing up, we plot Vi{N + |) versus Vs in Fig. 1, using Eqs. (7-11). The plot is the 
main result of this work. We can clearly see the single line of locahzation-delocalization 
transition, which has a relatively steep slope. The spacings between the adjacent lines 
are graduaUy decreasing. The slopes of lines are gradually increasing. The unexpected 
relatively large spacings appearing periodically is attributed by our limitation of the single 
level approximation. The spacings between adjacent lines, the slope of the special single 
steep line, the change of the slopes of small-slope hues, and the slope of the first line are 
controlled by the values of a, /3, 7, S, respectively. The starting points of the steep fine 
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and the first fine are determined by Sq and eV^^, respectively. Our plot of Fig. 1 shows 
the regular spacings between small-slope lines even for low electron densities. This is not in 
agreement with experimental results. It seems that the extended state Hamiltonian i^ext 
only valid for high electron densities. 

We can notice the anticrossings between the steep line and the small-slope lines. From 
the plot and Eq. (10), in order to study the anomaly of anticrossing, we should compare 
V^{N+l + l) in the region of e < D{N+1) withyt(-/V+|) in the region of e > D{N+1). We 
find that, if E{N + 1) + A(A^+ 1) - E{N) - X{N) is equal to or greater than E{N + 1) - E{N) , 
it is a normal anticrossing, and if otherwise, it is abnormal. In fact, we note A(iV + l) < X{N) 
for the cases of anomalous anticrossings. In the classical point of view, it is expected that the 
value of A(A^) is always increasing as becomes bigger. However, because of the quantum 
interference in the ground state of the system, it seems that the expectation value of A(A^-|-1) 
can be smaller than A(A^). 

In conclusion, we have considered a model Hamiltonian containing interactions between 
extended electrons and localized electrons. Coulomb blockade is found in energy calculation 
with only the extended state Hamiltonian. The general trends of the plot \\ versus Vs ob- 
tained by our theoretical study coincide with the experimental data. Including the localized 
state Hamiltonian, we explain the steep line observed in experiment. The possibility of the 
anomalous anticrossing is due to the interaction between extended and localized electrons. 
Our single level approximation introduced in this Letter is not applicable to the case of low 
electron densities. Perhaps, full calculation with more detailed Hamiltonian would be useful 
to explain the irregular spacings of the addition energies for low electron densities. 

The author is grateful to Professor C. K. Kim for introduction to quantum dots. 
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FIG. 1: The plot of the electron addition energy versus the side gate voltage is shown. Zoom-in to 
the part near = — 320meV is shown in the right. All quantities of energy dimension written in 
the right-top have the unit of meV. The dimensionless parameters a, /3, 7, 6 are chosen in order 
to obtain the similar feature of figure 1 (B) presented in Ref. 7. The values of A(n) related with 
C{Ncoie) is chosen simply as zero for all. The values of X{N) and t{N) from A?^ = 36 to 55 are 
arbitrarily given. A typical anomalous anticrossing takes place with A(39) < A(38). 



